High Tc superconductivity in MgB2 by nonadiabatic pairing 
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The evidence for the key role of the a bands in the electronic properties of MgB2 points to the 
possibility of nonadiabatic effects in the superconductivity of these materials. These are governed 
by the small value of the Fermi energy due to the vicinity of the hole doping level to the top of the 
a bands. We show that the nonadiabatic theory leads to a coherent interpretation of Tc = 39 K 
and the boron isotope coefficient qb = 0.30 without invoking very large couplings and it naturally 
explains the role of the disorder on Tc- It also leads to various specific predictions for the properties 
of MgB2 and for the material optimization of these type of compounds. 
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The field of high-Tc superconductivity is living an ex- 
citing time iQ]. New techniques provide in fact the 
possibiUty to explore physical regimes that were previ- 
ously inaccessible and superconducting materials which 
were often regarded as "conventional" BCS ones, as the 
fuUerenes, have proven to be real high-Tc compounds 
In this context the magnesium diboride MgB2, which was 
recently found to be superconductor with Tc = 39 K , 
is a promising material. The question is to assess whether 
MgB2 is one of the best optimized BCS materials or its 
superconducting properties stem from a novel mechanism 
of pairing and can be further improved in MgB2 or in re- 
lated compounds. In this Letter we would like to discuss 
some theoretical and experimental evidences that in our 
opinion point towards an unconventional type for the su- 
perconductivity, which we identify with the nonadiabatic 
framework. 

MgB2 is often regarded in literature as a conventional 
BCS-like superconductors, whose properties could be well 
described by the standard Migdal-Eliashberg (ME) the- 
ory. The high value of Tc is thus ascribed to the high 
frequency B-B phonon modes in the presence of a in- 
termediate or strong electron-phonon (el-ph) coupling A. 
LDA calculations find in fact A ~ 0.7 — 0.9 which, all 
together with the a representative phononic energy scale 
LUph ~ 650 — 850 K, is in principle able to account for the 
large value of Tc in MgB2 . However this picture is 
shaken by a series of facts. First, recent reflectance data 
are not consistent with a value of A strong enough to give 
Tc = 39 K HJl^. Second, the experimental determina- 
tion of the total isotope effect on Tc reported a boron 
isotope coefficient ae — 0.30 and a negligible magnesium 
isotope effect Preliminary indications suggest that 

this value of a cannot be explained by the LDA estimates 
of A ~ 0.7 — 0.9, but requires a much larger coupling 
A ~ 1.4 ||T^. We have solved numerically the Eliashberg 
equations to reproduce the experimental value of the iso- 
tope coefficient. We consider a rectangular Eliashberg 
function [a'^F{uj) = const, for 650 K < w < 850 K] as 
well as a simple Einstein spectrum with frequency luq- 
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FIG. 1. Critical temperature Tc as function of A for fixed 
value of boron isotope effect as = 0.30. Solid line corresponds 
to the rectangular Eliashberg function, grey region represents 
the solutions spanned by the Einstein model with frequency 
650 K < ojQ < 850 K. The dashed line marks the value T^ = 39 
K. 



The limiting values ojq = 650 K and ujq = 850 K of the 
Einstein model can be thus considered respectively as 
lower and upper bounds of a realistic Eliashberg func- 
tion. In Fig. |l| we show the critical temperature Tc as 
function of A for fixed value of a = 0.30. The correspond- 
ing needed Coulomb pseudopotential varies in the range 
fi* ~ 0.28 — 0.30 and does not depend on the specific 
Eliashberg function. We can see that a quite strong el- 
ph coupling is required to reproduce both Tc = 39 K and 
a = 0.30 with A ranging from 1.4 to 1.7. These values of 
A and fi* are thus even larger than the estimations of Ref. 
]lT[ , confirming and reinforcing the discrepancy between 
LDA results and the ME analysis of the experimental 
data (see also Ref. Q). Note that, contrary to cuprates 
and fuUerenes, electronic correlation is not expected to 
play a significant role in MgB2, and LDA calculations 
should be considered quantitatively reliable. 

This analysis therefore points towards a more complex 
framework to understand superconductivity in MgB2. 
An useful insight, in our opinion, comes from a com- 
parison of the electronic structure of MgB2 and graphite. 
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FIG. 2. Schematic band structure of graphite (top panel) 
and MgB2 (bottom panel). Grey region in top panel indi- 
cates the doping region achieved by chemical intercalation of 
graphite (±1 eV). 



These two compounds are indeed structurally and elec- 
tronically very similar. A main difference is the relative 
position of the a and tt bands with respect to the chem- 
ical potential /i. In undoped graphite the Fermi energy 
cuts the TT bands just at the K point, where the density 
of states (DOS) vanishes. Doping graphite with donors 
or acceptors, however, shifts the chemical potential /i of 
~ ±1 eV providing metallic charges in the system and 
a finite DOS This situation, on the other hand, is 

naturally accounted in MgB2, where /i lies well below the 
TT-band crossing at the K point and even crosses the two 
a bands (see Fig. ||, where a pictorial sketch of the band 
structure is drawn). Note that in the conventional ME 
context the only electronic relevant parameter is just the 
DOS at the Fermi level N{Q). From this point of view 
the difference between the superconducting properties of 
MgB2 with Tc = 39 K and intercalated doped graphite 
with Tc up to 0.55 K at ambient pressure is hard to jus- 
tify since both the materials show similar A^(0). Such a 
comparison suggests that the origin of the high- Tc phase 
in MgB2 should be sought among the features which dif- 
ferentiate MgB2 from doped graphite. 

A similar impasse was encountered in the ME descrip- 
tion of superconductivity in fuUerenes, which also share 
many similarities with graphite. Even there, LDA esti- 
mates of the el-ph coupling A were insufficient to account 
for the high Tc and for the small isotope effect. Such 
a discrepancy has been explained in terms of opening 
of nonadiabatic channels which, under favourable con- 



ditions fulfilled in fuUerenes, can effectively enhance the 
superconducting pairing [l^ . A key role is played by the 
small Fermi energy that in fuUerenes is of the same 
order of the phonon frequency, violating the adiabatic 
assumption (wph Ep). In this situation Migdal's theo- 
rem on which conventional ME theory relies, breaks 
down. The proper inclusion of the nonadiabatic contri- 
butions follows the framework of Ref. and leads to a 
new set of equations for superconductivity : 
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Z(a;„)A(w„) = Tc Fa (c^n ,UJm,Qc) ^^^"^^ IJm , 



(1) 

(2) 



where rim = 2 aTctan{EF/[Z{uj„i) ^^m]}, Z{ujn) is the 
renormalization function and A(w„) is the supercon- 
ducting gap function in Matsubara frequencies. The 
breakdown of Migdal's theorem strongly affects the "on- 
diagonal" Tz and the "off-diagonal" Fa el-ph kernels 
which include now vertex and cross contributions Jl^ : 

Vzil^n^i^m^Qc) = \D{uJn - iOm)[^ + AP(w„ , , Qc)] , 



Fa (ti'ri , , Qc) = XD{u!n — U!m)[l + 2XP{u! 



where -D(w„ — cOm) is the phonon propagator and fj, 
the dynamically unscreened Coulomb repulsion, to be 
not confused with the chemical potential. The vertex 
and cross functions, P{u!n,(^m,Qc) and C{uJn,^m,Qc), 
represent an average of the nonadiabatic diagrams over 
the momentum space probed by the el-ph scattering, 
parametrized by the quantity Qc- 

In the nonadiabatic context outlined above, the role 
of the a bands in MgB2 acquires a new and interesting 
perspective. Indeed the Fermi energy of these bands E^ 
is also quite small, E^ ^ 0.4 — 0.6 eV leading to 
^^ph/Ep ^ 0.1 — 0.2. These values, together the sizable 
A ~ 1, point towards a similar size of the vertex cor- 
rections AwphZ-Ep 0.1 — 0.2 and nonadiabatic channels 
induced by the breakdown of Migdal's theorem can be 
therefore expected to be operative. In this situation it 
is clear that the use of conventional ME framework can 
lead to inconsistent results and a nonadiabatic approach 
is unavoidable. The scenario we propose is the following: 

• MgB2 can be described as a multiband system with 
two conventional ME bands tt (with large Ep > 3 eV) 
and two nonadiabatic bands a {Ep ~ 0.4 — 0.6 eV). 

• TT bands can be in good approximation can be consid- 
ered as conventional. They could possibly contribute to 
the dynamical screening of /i* and to the static screening 
(Thomas-Fermi like) of the long-range el-ph interaction. 
They can also lead to the opening of a smaller super- 
conducting gap in the tt bands which does not probed 
directly nonadiabatic effects. 



2 



q 

E 

3 

3" 



0.4 



0.3 



0.2 



0.1 



0.0 





(a) 




- 


1 '1 

1 
1 




1 

1 

\ 

V 








0.0 0.2 0.4 0.6 

Q=q/2kp 

FIG. 3. (a) Momentum structure of the vertex function for 
a parabolic 2D hole-like band. Different curves correspond to 
different hole-fillings shown in panel (b). (c) Estimate of the 
enhancement of Tc for the nonadiabatic vertex theory with 
respect to the ME one. 



• High- Tc superconductivity is mainly driven by a- 
band states. The peculiar feature of such bands is the 
smallness of the Fermi energy which induces new (nona- 
diabatic) channels of el-ph interactions. Origin of the 
high-Tc superconductivity is the effective enhancement 
of the superconducting pairing as long as vertex correc- 
tions result positive [P(u;„, w™, Qc) > 0]. 

As seen in the last item, an important element in this 
scenario is the overall sign of the nonadiabatic effects, 
which governs the enhancement or the suppression of Tc- 
In previous studies we showed that the ve rtex function 
P roughly obeys the simple relation 
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P > 
F < 



v-pq/uj ^ 1 
v-pq/io ^ 1 



(3) 



where w is a generic exchanged energy involved in the 
scattering of order of tiiph, and v-p is the Fermi veloc- 
ity. In fuUerene compounds, the strong electronic cor- 
relation favours the small q momentum el-ph coupling 
[vpq/uj ^ 1] probing therefore the positive part of 
the vertex function P. 

In MgB2 the situation is deeply different. In fact the 
nonadiabatic regime in MgB2 is related to the closeness 
of the Fermi level to the top of the 2D cr-bands, and 
the non trivial dependence of the momentum-frequency 
structure of P on the filling has thus to be taken into 
account ||2^. To this regard, Eq. (||) is very helpful 
to illustrate this point since, for parabolic hole bands, 
f F cx where ^ is the chemical potential with respect 
to the top of the band. As /i is made smaller, the positive 
region of the vertex function will be enlarged and will 
eventually cover the whole momentum space. Hence, in 
MgB2 the nonadiabatic vertex diagrams are intrinsically 
positive in the whole momentum space regardless any 
electronic correlation. 

In Fig. H we show the numerical calculation of the 



momentum structure (panel a) of the vertex function 
P{u>n,ijJrmQ) [Q = q/^kp) for different hole filling of 
2D parabolic hole- like band(panel b). In panel a, the 
exchanged energy ujn — has been set equal to wo/2, 
where loq is an Einstein phonon representing the charac- 
teristic phonon energy scale. The structure of the ver- 
tex function is strongly dependent on the position of the 
chemical potential. In particular for almost filled band 
systems, as MgB2, the vertex structure becomes shape- 
less and positive (solid lines). In such a situation the con- 
tribution of the nonadiabatic vertex function is positive 
in the whole momentum space, and nonadiabatic chan- 
nels are expected to enhance Tc regardless the amount 
of electronic correlation. This trend is shown in Fig. 3c 
where the enhancement of Tc due to nonadiabatic vertex 
corrections is reported. The calculation of Tc follows a 
procedure similar to the one employed in Ref. [l9|] , where 
the vertex and cross functions are replaced by their re- 
spective averages over the momentum transfer and by 
setting ujn—u!m = wq. Note how, as moves towards the 
top of the band (panel b), Tc gets significantly enhanced 
by the opening of nonadiabatic channels already for val- 
ues of A consistent with the LDA calculations. Similar 
results were reported within the infinite dimensions ap- 
proximation [ p2| . 

It should be noted that the almost 2D character is an 
important ingredient for having a substantial value of Tc 
because the density of states remains finite at the band 
edge A 3D parabolic hole doped band would in fact 
lead to a DOS proportional to which vanishes as fi 

goes to zero. Additional effects can moreover arise from 
an intrinsic momentum modulation of the el-ph interac- 
tion. Low values of hole doping would in fact enlarge the 
screening length leading to an el-ph interaction peaked 
at small momentum transfer. A similar argument was 
proposed for instance in relation to copper oxides |2^] 
and, in principle, it could explain the refiectance data in 
MgB2 [^,|l0| . Both the argumentations can of course hold 
true and coexist in MgB2, explaining the high- Tc super- 
conductivity in this material as effect of a nonadiabatic 
el-ph pairing. 

We would like to stress that, once a bands are ac- 
cepted to play a key role in the superconducting pair- 
ing of MgB2 , nonadiabatic effects are unvoidably present 
due to the smallness of their Fermi energy. The onset 
of nonadiabatic channels can thus provide a natural ex- 
planation for the inconsistency between the theoretical 
values of A calculated by LDA technique (A ~ 0.7 — 0.9) 
and the high value A ^ 1.4 needed to reproduce experi- 
mental data Tc = 39 K and a = 0.30. 

Signatures of a nonadiabatic interaction can be found 
however in other anomalous properties of MgB2. The 
analysis of these features can provide further indipen- 
dent evidences for the nonadiabatic pairing and suggest 
precise experimental tests. 

Impurities and chemical doping. A remarkable 
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reduction of upon radiation-induced disorder has re- 
cently been observed in MgB2 , in contrast with An- 
derson's theorem. This kind of reduction in a s-wave 
superconductor has been shown to be one of the charac- 
teristic feature of a nonadiabatic pairing , as seen for 
instance in fullerenes ||2^. The experimentally observed 
reduction of T,. can be therefore a further evidence of 
nonadiabatic superconductivity. Similar conclusions can 
be drawn by the analysis of the chemical substitutional 
doping in MgB2. In fact, both electron |^ and hole ||2|] 
doped materials show a lower Tc than the pure stoichio- 
metric MgB2. It is clear however that the contemporary 
suppression of Tc upon electron or hole doping can not be 
understood in terms of band filling. We suggest a much 
more plausible scenario, namely that the stoichiometric 
disorder induced by chemical substitution to be mainly 
responsible for the reduction of Tc, with band filling as 
a secondary effect. Again, since nonmagnetic ion substi- 
tution does not break time reversal symmetry and An- 
derson's theorem in ME theory, a nonadiabatic pairing 
appears as a natural explanation. To test this picture 
the comparison with some completely substituted com- 
pounds would be interesting. 

Isotope effects. The detection of isotope effects on 
various quantities receives a crucial importance in the 
nonadiabatic framework since it directly probes the nona- 
diabatic nature of el-ph interaction. In particular it has 
been shown that nonadiabatic effects give rise to a finite 
isotope effect on quantities which in conventional ME 
theory are expected to not show it, for instance the ef- 
fective electron mass m* and the spin susceptibility 
X ||3^ ■ The actual discovery of an anomalous isotope ef- 
fects on these or other quantities represents therefore a 
precise prediction of the nonadiabatic theory which could 
be experimentally checked. 

New high-Tc materials. Interesting suggestions can 
come from the proposed nonadiabatic scenario in regard 
to material engineering and superconductivity optimiza- 
tion. According the analysis above discussed, a crucial 
difference between low- Tc doped graphite and high- Tc 
MgB2 is the upward shift of the a bands and their conse- 
quent cutting of the Fermi level. The study of the relative 
position of the a bands with respect of the tt bands, and 
of both of them with respect to the chemical potential 
appears therefore extremely interesting. In particular 
we would suggest that high- Tc superconductivity could 
be achieved in MgB2-like materials when Fermi level is 
lower but very close to the top of the a bands. On the 
contrary we predict no high- Tc superconductivity in the 
same family if compounds when i) Fermi level does not 
cross the a bands, ii) or where the Fermi level is very 
distant from the top of the a bands {E-p > 1 eV) and the 
system looses its nonadiabatic nature. Theoretical calcu- 
lations which can stimulate material engineering in this 
sense are in progress. A potential candidate would be the 
hole doped graphite as long as Fermi level could be low- 



ered to cut the underneath a bands or the a bands arisen 
by electrostatic effects. High level of chemical doping by 
acceptor intercalation was for long time unsuccessful in 
graphite as well as in Ceo since such compounds resulted 
unstable ]l3[ | . The recent discoveries of superconductiv- 
ity at Tc = 35 K in graphite-sulphur compounds J3^] and 
at Tc = 117 K in FET hole-doped fullerenes [g could 
thus both arise from the unifying framework of the nona- 
diabatic superconductivity. We thus encourage renewed 
work along these lines. 
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